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TECHNIQUE AND COMPUTER PROGRAM 
FOR CALCULATING PHOTOGRAPHIC FILM DENSITY VARIATIONS 
by Craig W. Qhlhorst 
Langley Research Center 

SUMMARY 


A technique and computer program have been developed to calculate the 
film density difference between the unexposed film border and any point on a* 
photograph. The program calculates the angle off the principal axis of the 
camera lens so that a correction can be made for vignetting and atmospheric 
backscattering. The program also plots the computed values as a function of 
position on the photograph so that a three-dimensional . picture is produced. 

Ranges of density difference can be predetermined, and the program will place 
each data point into its corresponding range so that the percentage of points 
in each rang§ can be calculated. 


SYMBOLS 


AA proportionality constant, mm/BCD 


AVERAGE number in BCD format representing the density of the unexposed 
film area 

BB proportionality constant representing the gradient of the linear 

density wedge, density /mm 

* 

BCD microdensitometer recorded film density 


CT 


represents the number of data points recorded from the border 
line trace 


Sum 


represents the summation of BCD numbers for the data points recorded 
from the border line trace 


XMM: 

ZZ 

abcd 

* 


calculated film density difference, density 

correction factor for vignetting and atmospheric backscattering 

difference in BCD numbers between any two, points 

angle Off the principal axis of the camera lens, degrees 



vertical component of <p, degrees 
c}>p.^ horizontal component of <J>, degrees 

0 camera horizontal half angle, degrees 

INTRODUCTION 


It takes considerable time and money to collect and analyze water samples 
to determine the health state of a body of water. Development of an instru- 
ment that could remotely determine the water quality of a body of water 
would be of great benefit to the effort of monitoring the water quality of 
this nation's rivers and lakes. 

The NASA Langley Research Center, in conjunction with other federal • 
agencies, has been working on a program to develop and evaluate instruments 
that will have the capacity to remotely determine various water quality 
parameters. One of the instruments being evaluated is a wide-angle lens, 
aerial photography system with appropriate lens filter combinations. The 
wide-angle lens is used to get both high spatial resolution and a synoptic 
view at the same* time. Film density of the photographs are then analyzed to 
see if there is a correlation between radiance intensity and certain water 
quality parameters. It has been shewn in reference 1, for example, that the 
film density difference between the unexposed film border and a point on a 
photograph can be correlated to chlorophyll ji concentrations. 

Data obtained using the wide-angle lens require a mathematical correction 
to the film densities to compensate for atmospheric backscattering and lens 
vignetting. The correction needed is a function of location on the photograph 
(angle off the principal axis of the lens). The principal axis is assumed here 
to be the same as the perpendicular vertical intersecting the center point of 
the photograph. Each data point will have 'a unique angle, so a method to 
calculate this angle is needed. A microdensitometer, an isodensitracer , and 
a magnetic tape encoder were combined to digitalize the film densities. This 
paper presents a technique and computer program to convert digitalized film 
densities into a number representing the film density difference between the 
unexposed film border and a point on the photograph. 


PROBLEM TASK DESCRIPTION 


A National Instruments Laboratories' Isodensitracer and Magnetic Tape 
Digital Encoder System are used to digitalize film densities obtained 
from a Joyce, Label and Company Microdensitometer. These are recorded in 
Binary Coded Decimal (BCD) format. The problem considered in this paper is 
that of calculating the film density difference between the unexposed film 
border and any recorded data point. 
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TECHNIQUE FOR COMPUTING FILM DENSITY DIFFERENCES 


The Joyce, Lobel and Company Microdensitometer , Model Mark III C, has 
been modified electronically to utilize the automatic interval stepper of the 
isodensitracer so that the microdensitometer can automatically scan a full 
photograph* The microdensitometer uses a double-beam system in which two 
light beams from a single source are alternately impinged on a photomultiplier 
tube* The photomultiplier alternately sees the beam going through the speci- 
men and through a reference linear density wedge* thus providing an electrical 
signal which causes movements of a servomotor to bring the linear wedge into 
optical balance with the specimen beam* A pen attached to the linear wedge 
carriage plots the position of the wedge resulting in a plot of the density. 

The vertical distance that the tracing pen moves, (see figure l) , equals the 
distance that the linear wedge is shifted to optically balance the reference 
and specimen beam. Multiplication of this vertical pen distance by the 
slope (BB, density /mm) of the linear density wedge will give a film density 
difference value. 

The Magnetic Tape Encoder System is mounted on the rear of the 
microdensitometer and attached to one of the pulleys driving the linear 
density wedge. A multi commutator is mounted on the x-axis of the flat bed 
recorder. This commutator provides 5, 10, 20, 50, or 100 contacts per revolu- 
tion and is so mounted to provide approximately 10 revolutions for one traverse 
of the recording table on the X-axis. At each point of contact, density datum is^ 
recorded. The encoder provides a BCD output corresponding to one of 175 posi- 
tions evenly spread out along the length of the linear wedge; thus, there is a 
direct proportionality between the magnitude of the difference in BCD numbers 
between any two points and the vertical pen distance between those points* 

Vertical Pen Distance = (AA)(ABCD) 

The proportionality constant (AA) has been hand calculated and stored 
in the program* (See Appendix A for the calculation of AA.) A mathematical 
computation subtracting the recorded BCD number for a specific data point 
from the BCD number representing the unexposed film border (AVERAGE) will 
give the ABCD number for that data point* AVERAGE is calculated by summing 
up the BCD numbers of the border trace and dividing by the total number of 
points in that trace , 

AVERAGE - Sum/CT 

As previously stated, multiplication of the vertical pen distance by 
constant BB would give the film density difference for a specific data point. 

Use of a wide-angle lens requires that a correction for vignetting and 
atmospheric backs cattering be made* The correction must be applied to the 
vertical pen distance before multiplication by BB. The correction factor (ZZ), 
can be calculated and applies to each data point, provided the angle <j) for 
each point is known, (see fig* 2)* For the photograph shown in this paper, 

ZZ was set equal to cos^cf)* The angle can.be broken up into its horizontal 
and vertical components with 
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2 2 2 
tan <p = tan ' <j^ + tan 

The angles <b . and d>. . can be calculated since the number of points 
r pi Aj 


recorded per scan and the number of scan lines are known. 

Camera horizontal half angle 
^pi “ l/2(number of points per scan) 


X p. 


diA • ” Camera horizontal half angle 

j l/2(total number of scan lines) j 


where 


= ith point from the center of the film 


A. 

J 


= jth line from 


These angles are calculated by the 
the user, an'3 stored by Data Statements 


the appropriate * tan <(> 

together, and calculates 
density difference are: 


and 

4 >- 


tan^<t> 

Y Aj 

Thus, the 


the centerline 

2 

user, converted to tan values by 
in the program. The program calls 

values for each data point, adds them 

steps for calculating the film 


1. Calculation of AVERAGE 

2. AVERAGE - BCD = ABCD 

3. (ABCD) X (AA) - Vertical Pen Distance 
Calculation of (f> 

5« Calculation of ZZ 

6. (Vertical Pen Distance )/( ZZ) = Corrected Vertical Pen Distance 
7- (Corrected Vertical Pen Distance) X (BB) = XMM 

L 


PROGRAM DESCRIPTION 


Program CONVERT is written in Fortran IV for the Control Data 6000 series 
computers. It uses Langley Research Center Library subroutines PSEUDO, 

CALPLT , AXIS, LINPLT , N FRAME; and Library functions SQRT , ATAN, and COS. 

The program also uses Subroutine ENCOD, which is not in the Langley Computer 
Library but was developed by Langley personnel for use with the Magnetic 
Tape Encoder. 


Program CONVERT 

Program CONVERT is used for the input and output of data. It initializes 
the variables. . Subroutine ENCOD takes the information off the tape, and the 
program calculates the corrected film density difference values. Each value 
is placed in a range of density differences by the program so that the 
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percentage of points in each, range can be calculated. The user specifies 
the desired density difference ranges. Program CONVERT also plots the 
computed values as a function of position on the photograph so that a three- 
dimensional picture is produced with the third dimension being the density 
difference . 

A specific order of data input onto the recording tape must be followed 
for the program to work properly. The specific order needed is explained in 
Appendix B, and a program listing is given in Appendix C. The flow diagram 
of program CONVERT is shown on pages 6 and 7« 


PROGRAM USAGE 


The program is run on the Control Data 6000 series computer under the 
Scope 3.0 operating system and requires a field length of 60 ,000g storage . 
locations. The Central Processing Unit (CPU) time will depend on the number 
of data points recorded. A CPU time of 150 seconds will usually be enough 
to analyze one 70 mm size photograph with 50,000 data points recorded. 


Input Description 

Most of the input data comes from the user submitted tape. Other input 
data terms are as follows; 


PPP 


CFX 


CFY 

AA 

BB 


array of Kodak step wedge densities corresponding to the wedge 
steps that were traced over by the microdensitometer, inputed by 
a Data Statement 


array of the tan <j> 


Pi 


values for each data point in a line, 


inputed by a Data Statement 
2 

array of the tan $ values for each line, inputed by a Data 
Statement ^ 


constant, m m/BCD 
constant, density /mm 


Desired density difference ranges — entered by IF Statements. 

Output Description 

The output for program CONVERT consists of these elements: 

1. Five lines of machine setting and other non-density data. See Table I 
for sample listing. The settings are recorded by the Fixed Data Switches on 
the Magnetic Tape Encoder (see Appendix D). 
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continued Flow Diagram 
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2. Table of step wedge density versus assigned BCD number. Analysis 

of the BCD numbers will determine if the microdensitometer is working properly. 
There should be a constant difference, -1, between any two BCD numbers. (See 
Table II. 

3 . A counter (y-number) specifying what trace the upcoming data repre- 
sents. A zero means that the border trace is printed next. See top left 
corner of Table III. 

i+. The border trace is printed in BCD format so that the user can 
check the density variation. From these data, AVERAGE is calculated and 
printed out. See Table III. 

5 . Calculated density difference values for each scan. There are 
twelve data points printed per line of printout. See Table IV. 

6 . Table of density difference range, number of points in each range, 
and the percentage of total number of data points in each range. See 
Table V. 

CC represents total number of data points 

EA, EB, EC,* ED, EE, EF, EG, EH represent the number of data points 

in specific difference ranges. 

7. Computer plot The plot is a three-dimensional picture of 

the calculated film density differences. 

Each line is plotted from right to left. The 
actual scales allow for each point on a line 
to be plotted 2.|>4 mm ( 0.1 in.) apart, and 
the starting point of each successive line is 
shifted 0.762 mm (0.03 in.) to the left. See 

figure 3. 

% 

Machine Setting Terms are as follows: 

A height of light spot on the aperture plate 

of the microdensitometer , recorded as two 
BCD numbers 

B width of light spot on the aperture plate of 

the microdensitometer, recorded as three BCD 
digits 

G width of li^it spot on the photograph, 

recorded as three BCD digits 

condenser focuses the light beam on the photograph, 

number represents the lens power, recorded 
as two BCD digits 



wedge 


commutator 


dif* control 


pen damping 


write-out arm ratio 


pen.^stop multiplier 


spec, step multiplier 


objective 


magnification 


resolution 


film number 


cal. strip number 


represents density- range of the linear wedge, 
recorded as two BCD digits 

represents the number of segments that the 
linear wedge is broken up into by the 
isodensitracer , recorded as two BCD digits 

an electronic control which varies the amount 
of feedback signal to the table drive unit, 
recorded as two BCD digits 

controls the speed of the tracing pen when 
the microdensitometer is used in manual mode, 
recorded as two BCD digits 

ratio of recording table movement to 
specimen table movement, recorded as three 
BCD digits 

number of 125 micron-steps that the iso- 
densitracer recording pen takes between 
successive scans, used when the system is 
set to run in automatic mode, recorded as 
two BCD digits 

number of 25 mi cron- steps that the specimen 
table is moved between successive scans, 
used when the system is set to automatic 
mode, recorded as two BCD digits 

focuses light from the specimen table onto 
the aperture plate of the microdensitometer, 
number represents the lens power, recorded 
as two BCD digits 

built into the microdensitometer, number 
represents lens power recorded as two BCD 
digits 

number of data points recorded by the 
magnetic tape encoder system per inch of 
recording table movement, recorded as 
three BCD digits 

number assigned to distinguish the photograph 
from others , recorded as five BCD digits 

number specifying the Kodak step wedge used, 
recorded as five BCD digits 
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step low density- 


number of the lowest Kodak density step 
traced over, recorded as two BCD digits 

step high density number of the highest Kodak density step 

traced over, recorded as two BCD digits 


SAMPLE CASE 


"A positive transparency of a 70 mm photograph taken above Maryland Point 
on the Potomac River in October 1972 was used for the sample case. (See 
figure 1.) The picture was taken by a Hasselblad 500 EL camera at an altitude 
of 3,01+8 m (10,000 ft). The camera had a 1+0 mm focal length and a Wratten 
89B filter was attached. An F-stop of 8 was used with a shutter speed of 
1/250 seconds. The positive transparency was overdeveloped to bring out the 
features in the water. The white area covering the top third of the photo- 
graph is land. The rest of the picture is water. Ground truth data have shown 
the water to be heavily concentrated with blue-green algae. Analysis has shown 
that the white sections in the water are areas with sufficient algae density 
to produce chlorophyll concentrations greater than 3*+ UgA • 

Density data on the photograph were measured by the microdensitometer 
from left to right in the horizontal direction and from bottom to top in the 
vertical direction. The system was set up so that 1+06 data points were 
recorded from each scan with 101 scans evenly spaced and covering the whole 
photograph . 


Sample Input Data 


PPP 0,05, 0.20, 0.36, 0.50, 0.6I+, 0.79s 0.94, 1.09, 1.23, 1.38, 1.53 

0 

CFX the array of tan $ ^ values inputed as printed out in the program 

listing (Appendix C) 

2 

CFY the array of tan <J> values inputed as printed out in the program 

A J 

listing (Appendix C) 


AA = 1.1530 mm/BCD 

3 

BB = 0.0130 density change /mm 

Density difference ranges <0.11+3, 0.11+3 - 0.195, 0.195 - 0.260, 0.260 - 0.299* 

0.299 - 0,51+6, 0.51+6 - 0.650, 0.650 - 1.1+00, 1.1+00 


Sample Output 

The printed 'output is shown in Tables I through Table V with the 
computer plot shown in figure 3. 
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CONCLUDING REMARKS 


A technique and computer program have been developed to calculate the 
film density difference between the unexposed film border and any point on a 
photograph . The program calculates the angle off the principal axis of the 
camera lens so that a correction can be made for vignetting and atmospheric 
backscattering. The program also plots the computed values as a function of 
position on the photograph so that a three-dimensional picture is produced. 
Ranges of density difference can be predetermined, and the program will place 
each data point into its corresponding range so that the percentage of points 
in each range can be calculated. 
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APPENDIX A 


CALCULATION OF CONSTANT AA 


Plotting of a density profile of a Kodak Step Wedge (see figure 4) 

■while at the same time recording the data on tape enables the constant AA 
to be calculated. The constant is computed by measuring the vertical distance 
between any two steps of the wedge on the profile and dividing by the differ- 
ence between the corresponding BCD numbers. The BCD numbers are taken from 
Table II . 

Example : 

The distance between steps 3 and 7 of figure A1 = 62.2 mm 

The corresponding BCD numbers are 97 for step 3 and 151 for step 7 • 

ABCD = 151-97 = 5b 


62.2 mm 1.1 518 mm 

5 Vci> " iECD 

The value of 1.1530 mm/BCD assigned to AA in the sample case came about 
from the averaging of ten such calculations using the F-lh0» 2.b density wedge 
but should be good for any linear wedge used. 
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APPENDIX B 


ORDER OF DATA INPUT NEEDED FOR PROGRAM CONVERT 


To insure that program CONVERT works properly, a specific order of data 
input onto the recording magnetic tape should he followed* 

First - Five lines of machine setting data should he inputed onto the 
tape* Use the fixed data switches* 

2nd - A trace of a Kodak step wedge should be inputed next 

3rd - A border trace should then be inputed 

After these three steps have been followed, the program is now ready' for 
these film data. 

When the user is finished with one photograph, a set of fixed data should 
be inputed onto the tape with the first switch being set to 3 and then another 
set of fixed dat^L placed in the tape with the first fixed data switch set at 9* 
See Appendix D for an explanation of this first switch. When digitalizing a 
second photograph on the recording tape, the first and second steps from above 
can be left out. 

When the user is finished putting data on a tape, an END OF FILE marker 
(EOF) should be placed on the tape. The EOF signifies to the computer ‘that 
there is no more data on this specific tape and thus terminates the program. 
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APPENDIX C 


Program Listing 



EXIT. 

ONLOAD (TAPE12) 


I* LTiilfl'M LQ *11 


PEN j$ A~TW0~0IHEN51QNAL ARRATt CONTAINS DECODED DENSITY INFORMATION 

TWn»PtMFNSIQNAL ARRAT* CONTAINS DECOPEP F 1 X ED DAT A INFORMATION 

CPA AR&AY OF DATA PT5. f AxFN FTOH A HICRa _ DENSJTQMETER TRAtt I 


OF A KODAK STEP WEDGE 


D ONE-OI MENS ION ARRAY CONTAINS SAME INFORMATION AS IDEM 



ill iM*i Mnw-wnifMJ 


XKK ARRAY the LOCATIONS ON THE HORIZONTAL AXIS OF THE GRAPH 


PAPER WHERE THE OAT A POINTS ARE TO BE PLOTTED 




1 RO ( 51 2 ) 9 XKK (512) * XHM <512) *CFX <430 ) t CFY ( 1 1 0 ) 



tograph 


DATA (CFX <n*I*lfl25>/0«4727l* 


•44917 t 


_ 

LMmMSm5i5SBa^EaaEEmiiEaBa^^^ ^^^^ ^BI 


30 . 35933 » 0 . 35A5 1 1 0 . 3A975 . 0 * 3450A * 0*3A037t 






Ul.tl 


50. 28796. 0. 2839 1. 0, 2?987 ♦ 0. 27587, 0.27 192 « 0,267997 b. 26412.0, 26029. 

60. 256*9.0. 25274,0.24901 .0.24S 34, 0.241 70.0.2381 0,0.23453.0.231 00, 
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TO. 22751 ,0.22405, 0.22062, 0.21724, 0.21390, 0.21 058, 0.20731. 0.20406, 
80. 20085,0. 19767, 0.194S3, 0.19 142, 0.18834, 0 .1 8529 ,0. 182gg * 0 iil311i_ 

on 1 17145*0. 17054*0.16770*0* 16 486*0.16207 g _0 *1593. 1 1 0 * 

~l^t~1^3fi9^0^15T2?«^T4H56^0« 14595* o! 14338,0 1 1408 3^0 .J_3B3J L jJL»A3iiX^i- 


80. 13333 ,0.13089,0. 12848,0. 12610,0. 123 72,0. 1 2140,0. 1 1909. 0. 1 1662j _ 
rn. ii 4 «? 7 .n-i 1233 . 0 . II Q 14 . ~fl. i»796. 0. 10582, 0 . 10 3 70 » ( L. l til 60,0 . 0995 2.. 


rn. i i467.n. n?33.n. n oi . i n r*t>. u. * uaarm . i iiji n»« » » 

DO*OQ72HTnt»954ft!nIn4347.oIo914 9. 0.089S4.0.087 63j»jjuOi^7-3x P*^^^^^* 

FO.na2no.o.OBOl7.O.O7B37.O.076 6 0 , 0.07484,q.Q73 I O,0.Q7t39,0.P697 1,- 

ro. 0^804. 0. 0664 1) , P.0647fltfl . 06 j l B ^. 0ftljSl,0 . Q60g6j> 0 . 05^5-1,0 . 0570?/ 

~ J pftifA(CrXlT)«l = 126, 250 >/OT05553 ,QTo 540 7 ,0.05262, Oj 051 20 , 0.04^ 8g^ — 

— — : 1 ■ — _ _ ; I. . . . A a < -i i i n a l o . a niiACll.n.n IS , 


1 0 .04842,0.04706,0.04572,0.04441,0.04311 ,0. 04 1 B^jO gOMgg j 0*03935 . , 

?n . (1381 4.0.03694. 0 . 03577 , 0 . 03462.0 .03349, 0 .03237 ,0 »Q312B_, 9 * 

3ft:°?9u:o:^an:o:0^nr6. 02612^. 0^^0.02419, 0. - 0 2325,0.022^ 


'4n.A2l44 .^. 02056. 0.0 1970.0.01886,0. 01804, 0 . 01 72 4,0 .01 6 46,0.01569, 

<^n.n 140^.0.0 1422. 0 .0 1351 ,0.01282,0.01 2 1 5,0.01150,0.01^86,0.01025, 

60.00965,0.009^7,0.00851 , 0. 00797 ,0.00 744,0, _00693^0^00645»0.0q|98»_ 


-"TO. 00552, 0.0 0509, "0.00467, 0.00427, 0 .00389,0.00353,0 . 0031 9, 0.002.8^ 

80.00255.0 . 00226,_9 . 00198 , 0 < 


8n.(l0255 .0.00226.0.00198,0.00173,0^ o oi4v.tf.uuicrt»tyv»v'T».yvv^ .4. 

4o!o6&7i I o!oo05 6lo!oo643. 0.00032.0.00022,0.00^14,0.0000^ ^^00^ 

AO.OOdO) tO.OOOTO.O. 00001 ,0.0000470.00008*0. 00014, 0^^0^2^0_^00032»^ 
^QToO O^TO^OOSt^. 00071.^ 00 088. 0.-00107,0.00127, 


BO. 00043 ■Q*OOOb b*0«UUU f i fU* U UUPQ* r t **»**- -** 

~7!g700 196. 0. 00226, 0. 00255, 0.00286,0.0031 9,0. 00353, 0.00389j0^00422j_ 


D 0 . 0 0467 >0. 00509^ 0. 00552, O.QQ59B7o. 00645,0. 00 6 93, 0.00744,0. 00 79 7j_ 

To. no 851 . o. 0040 7.o. 00965. 0.0102 5,0 . 01086, Q . Q l l 50 id. 0 1 21 5 ,0 . Q 1 2S 2» 

Folo 1351, olo 1422,0 io 1495,0. 01569.0, 01646,0.01 724,0j 0l804. 0 j0l 886/ ^ 


Hat A ( CFX < 1 ) . I *25 1 . 375 ) /0 ,1Tl 97 0 , 0 ,020 5 6 ,0,021 44 . 0 , 02^34, 0 , 0 2 31^3 




3Fio4T84. 0.04311. 0.04441, 0.04572. 0.04706. 0.04842 , 0.0498 0 ,0. 051 2 gj_ 


40. 05262. Q.Q5407. 0.05553. n.05702,0. 0 5853, Q »(j6006tQ,0 6 l 6 1 1 0,063 1- 9 - l 
50. 06478, o! 06640, 0*06804, 0.06971 >0.07139,0.0731 0,0_»j 


.074B4, 0.07660, 


60.07837,0 .080 17. 0,08200.0. 08386. 0. 08573. <K 08763 ,0.0 895_4. 0.0 91 


"7ol 09347. 0.09546.0.09748 , 0.09953, 0.10160,0, 10370,0, 1058g,0,lQ7*6j _ 

an! i mi4.nl 11 233.0 I 11457. 0.11 682, 0.1 1909, 0 . 12140,0. 12372, 0 . 12610, 

^rrr28^rrT3d84;o. 13333.0. 13579, 0. 13831. o.iA0-83>QjA»m»o J ^^ 




t 164BB,p, 16770^. 
.18834,0.19142, 


80.1705470 


CO. 19453. 0 .19767 ,0.20085.0 .20406.0.20731,0.21 058, 0 .21 390 ,_0 « 21724. 


22405.0,22751 ,0 . 23100.0.23453,0.238 1 0,0 


DO. 22062.0 JC 

Ed. 24yol, 0.25274, 0^25649, 6 
FO. 27987, 0.2839 1, 0.28798,0 


.24170,0.24534, 

.2Tm70T2758T7 


726029, 0 .26412.0, _ 

7r92TT7o7?^527. 0.30048, 0 . 3047 4, 0,30964/ 


?0 -37405. 0.37907. 0.3841 3.0 


- j3892St0t39 44g t 0 - 


.35933 .0. 

.39964,0. 


40492,0.41025, 


TAN**? VALUES FOR EACH LINE STARTING AT BOTTOM OF PHOTOGRAPH 


.0.43777,0.44344,0.44917,0 .45496. 




?n . ??H34. n.?l 452.0^20 125,0718852,0.1 7636 ,0.16472,0.15359,0.14294, 
30. 13274. 0.12303, 0.1 1376*0.1 0493, 0.09652*0. 08851, 0.08089, 0.07368, 
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40. 06685. 0. 06039. 0.05429. 0-04854, 0.04315, 0.0381 0,0. 03339. 0.02901, 

50.02494.0. 02120.0.61777. 0.0 1466.0.01 185. 0_,_0 0?_35.0. 00715, Q._00jg4_i 

-6 070036^, 0.00 233,0. 00 1 31^0058^0 . 00015, 0^-^* 




JL 

_c_ 

c 


777 .n.ft 71 ? 0 .(l . 02494 , 0 . 0 2901,0 .03339, 0 . 03 B 1 0,0. 043 15,0 .0485 , 

90. 05429, o - Q4Q -39. 0.06685,0. 0736 8, O.Oft089 T O t pB85ltO.P^5g, , 0, | 0493, 

- Ir?l^:“: - l^ n 3loll32/4 ,0.14294 ;o:i5359;o. 164/^.0.17 63 6. 0.18852. __ 

80.20125,0.21452, 0 . ??8 36. n . 24279. 0.257B6.0 .27362, 0.29 0 00, 0. 3Q7Q9. 

ro. 32490. 0.34348. 0.36290, 0.3 8 309, 0,40414, 0.42610 ^. 0, 44893, 0.47£ ZL. 

09*0.00000/ _ 

<Q 1 IN^TUM Y SET EQUAL TO 7ER0 AND THEN REPRESENTS THf Qf 

1 TIMFS THAT SUBROUT INE ENCQO IS CALLEP 

X KK INITIALLY SET TO 2ER0 . REPRESENTS THE NUMBER Of UNf TRACES TA«EN- 
— ™ Afi l I rn^i FR TQ SE ^~l^^GiU*irTO^ER 0 J^HEM_^PPRQPR1A1£_ 


~C C = TOTAL NUH h ttH O F DlMIAUztU umiA POINTS M , p , ■ 

Ea.EB Etc. = NUH6£ff ~ 0P DlSlTALlO EP DATA POINT* IN VftRIQU* PEW U fF RANOfcS 


<0=0 
lff T~ i ^ 


TT 


mKO.EQ.Ol GO TO 1 1 0 

3) 


CALL CALPLf 14. .-3., 
110 NO=0 


CC=0. 


EH=irr 


HC=0~, 


£D=0. 


EE=0 . 


EF=0. 


£H=UT 


10 KO=KO, 1 


Hurt mumrTb OF RECORDED QAtA PQ I NT 6 IN A LINE 
CALL EnCODUCNT.O. lDtN.lQ,N6C8t5,NC0Ut) 


PRINT 25 


25 FORMAT <1H01 
IF mCo&L.EQ.Oi GO TO 46 


IF(NCoPE7E0.2T GO TO 30 
1CNT=NSCPTS 
NXPT=NSLPlS.r 


C CON VErS16 N~6F P TWQ-DIHENSIONAL array to a one DIME N SIONAL ARRAY 

— 601 om = iDtmwiT ~~ 


IF (KQ .NE.6) GO TO 26 . - CT — 

a 99 ? r.uurMT OF a urn NO. TO A SPECIFIC KODAK MEP HEDGE D ENSITY 


00 2 NA=S.NXPT,36 


<e=K8* 


CPA(KB>=U INAT 


2 CONTINUE 


PRINT 88 
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88 FORMAT (1H1 »«OENSlTY* 1 0K*BCD NO.«> 

PRINT 89« CPPP ( I ) «CPA CI)tI~l»KB) 

89 FORMAT (1H0»F6.2*12X*F3,0) 

GO TO 10 

26 IP=M 

NOsNO* 


IFCNO.EQ.l) KK=0 

PRINT 22«KK . 

22 FORMAT ( 1H0 »AX»*V=*»I10//> 

00 I 1*1 ♦ ICNT 

NEEDED TO NUMERICALLY ORDER THE BCD NUMBERS 


IF CD 1 1 1 .LT»CPA ( 1 ) ) P<I>»D(I)* 1 7 5. 

IFCKK.EQ.O) GO TQ I 

CALCULATION Of DELTA BCD NUMBERS ; 



m 

IlflM 


Nr.LT. 10) GO TO 350 

IF (KK,GE* 1 > GO TO 902 
PRINT 200.1DU)»Igt»ICNT) 


200 FORMAT <6X»12F10«2) 

CT = 0 


SUH = 0 



AVERAGE-5UM/CT 


PRINT 900 


PRINT 90 1 t AVERAGE 


901 FORMAT <52XfFJ0^> 



CC=CC*1- 


XMM ( I 


800 CONTINUE 


IF (KK, « L T » 102) GO TO SOi 


iezhu 


351 FORMAT THESE VAL UES HAVE NOT BEEN CORRECTED' 8 ! 

C OFF-CFNTER ANGLE AND CORRECTION FACTOR CALCULATIO N 

803 DO 602 1=1* ICNT 


ORIGINAL PAGE IS 
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CF=CFY(KK) 

CFF=CF*m 

IF (KK.GT .101 ) GO TO 330 

TA~NPHI=SQRHCF»CFF) 

PHI=ATAN(TANPHH 

ZZ*COS<PHU.»»£_ 

60 TO 331 

330 ZZ»1« __ 

CALCULATION OF VERTICAL PEN DISTANCE WITH CORRECTION FACTOR APPLlE 

331 XHM<I)=XHH(I)/Z2 

802 CONTINUE 

3 . -- 


20 3 FORMAT <58X«13H0ELTA DENSITY) 





DO 80 



201 roRMftT(6Xtl2riQ»<») 

IF(KK.GT.lOl) GO TO 10 

~C ASSIGNMENT OP EACH DATA PT* TO A DENSITY DIFFERENCE RANGE 

D0~3frl ~I = ltlCNT _ 

IF ( XMM (I) »G£ . 0 . 1430 ) G0TQ3Q2 

EA-EA4 1 


302 IF(XMMU) ,GT. 0.1950) GO TO 303 



303 IF(XHh(I> «GT. 0.26QQ) GO TO 304 


304 IF(XMM<I > .GT. 0.2990) GO TO 305 


ED=ED*l 


305 IF(XMM( I) .GT. 0.5460) GO TO 306 


EE=EE* 



301 


309 £H*EH+1. 

301 CONTINUE 


IN 

XKKU )=0. 


00 4 1-2* ICNT 


l 

4 CONTINUE 

XMM ( I CNT ♦ 1 ) =0. 
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220 CALL L1NPLT(XMM<)(KK«ICNT<1«0»0»0*OI - 

CALL CALPLT («10«»03»~3> - — 

_ _ I F (KK , EQ » 1 0 1 L . 60 TO 3Q7 

GO TO 10 

CAlC UlAT ION Of PERCENTAGES OF TOTAL PTS. TOR EACH PEN D IF RANGE 
EAAj ERA ETC. REPRESENT THE PERCENTAGE OF THE TOTAL NU NBEg Of PAlA_ 
i POINTS In SpEcIFiC DENSITY &IFFERENCE RANGED 


EBA= (EB/CC) *100 
ECA= (EC/CC> B 100 
EDA= <EO/CC) » 100 




EFA= (EF/CC1 »100 

EGA= (E6/CC1 * 100 

forhat H £ta$enents for QEnSIT Y oiffErEnce r ange table 

PRINT 310 

310 FOPMAT(/60X«*OEN DIF RANfiE»« IPX »*N0« PTS,», _ 


1 10X»«PERCENT AREA* > 

PRINT 3H«CC ■ 

311 FORMAT </20X.2HCC«63X oF6,,0> 

PRINT 3 1 2 » E At £AA 

312 FQRHAT(20X9 2HEA«*>0X. < >.LT.0.1A3»d*X.F6.0» l»X> 

PRINT 3139E89ETa _ _ZIZ - 




F6.2> 



DO 21 1=1 « 12 

IF(10(I > »EQ tt 131 iod>=o 
21 CONTINUE 


PUGinai; page is 
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CALL CAL PL t <0» dOo *999) 

STOP 

END 

5UbR0UtIn£ EnCoD I iCNt »k5kTp * DE'N* ID»NSCPTS*NC00E) 

INTEGER OEN(500> *ID<12> «A(10A9> 

INTEGER CON (15) *LM ( 5) «LL ( 5) 

HFDl =000000000000000000060 


MASK 2=000000000 0 0000000001 7 
MASK 3=077007700770077007700 
IE ( KSK I P »EQ. 0) GO TO 999 

DO 998 J=1«K5KIP 

950 ICNT = IVAR<12,A,5» 1050) 

IF (ICNt .EO. 01 GO 10 998 


GO TO 950 


998 


CONTINUE 
<SKIP = 0“ 
999 CONTINUE 


ICNT = I VAR (1 2 * A * 5 * 1 050 ) 
FORMAT t2X*MtNTlS*I10) 


T5F 


IF (ICNT *EQ» 0) GO TO 77 
A ( 2 > = A ( 2~) iO. MASK3 
DECODE (10 * 2 0 » A ( 2 1 ) 


1 LL (MM) »MM= 1 ♦ 5 ) 


20 FORMAT < 5 ( 1 X * R 1 ) ( 


LI =MFD1 iAi LL(3) 

IF < L 1 . EQo O.OR. Ll*EQ.AQ8) GO TO 70 

1=0 ' 


LfcQ 

1 DO 1 0 J=l*3 

I = I £l_ 

IF <1 .GT. ICNT) GO TO 14 

A< I ) =A( I ) .0. MASK3 

DECODE ( 10 *20 « A ( I) ) <LL<MM) *MM=l*5) 

DO 99 K«l*5 

M=S* < J»n *K 


~CON(M)sL L(K> .A. H ASK 2 



99 CONTINUE 


10 CONTINUE 

27 LJ=M/3^ 

IF (LJ «EQ« 0) GO TO 15 


00 11 J=1 «LJ 
K=3*(J-1> 


L C L ♦ 1 _ 

DEM (I > s 1 00* CON (K» 1 ) *10 »Cf)N ( K ♦ 2 ) ♦ CON < K ♦ 3 ) 
11 CONTINUE , 

IF <LJ .LT. 5)~G0 TO 15 

GO TO 1 

u j(mu6*i(m 

101 FORMAT 1//10HREAD ERROR ///) 

” 15 NSCPTS=L-1 - 
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10(L) =LL (K1 .A.~ MASK 2 

50 CONTINUE ' ’ 

' 60 NC0DE=2 

— 

T7 MCoor=o ' 

RETURN -7" 



APPENDIX D 


FIXED DATA SWITCHES 


The Fixed Data that is input ed onto the recording magnetic tape is 
recorded as 12 separate Binary Coded Digits . 

Switches 


1 

2 i 

3 

4 

l 5 

6 

7 

8 

9 

10 

11 

12 


BCD digits, 0-9, can be recorded by each of the twelve switches. Switches 
two through twelve can be used to record non-density film data, i.e. machine 
setting data, etc. 

Switch 1 is used as a code: 

If Switch 1 

ft 

equals 0 - means that the rest of the fixed data on that same line is 
non-density film data 

equals 1 - means that the next set of non-fixed data represents the 
microdensitometer trace over a Kodak Step Wedge 

equals 3 - tells the computer to go to Fortran Statement number 307 

which starts the computation of the number and percentages 
of data points in specific density ranges up to that point 

equals h - means that the following data are information collected from 
the photograph of interest 

equals 9 - signifies the end of recorded data for a particular photograph * 
tells the computer to go to Fortran Statement number 107 which 
resets all counter variables back to zero for the start of a 
new piece of film. 
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TABLE I.- MACHINE SETTING 


COPi 

A 

8 


G 


HON DENSE* 

0 

6 5 

1 5 

0 1 

5 



l 6 0 

* 

core 

WEDGE 

COMMUTATOR 

OIF CONTROL 

PEN 

DAMPING 

WRITE-OUT RATIO ARM 

0 

2 4 

4 6 

0 _ 5 

i Q_ 

t 

1 

J 

G 

O 

* i 

none 

PEN STEP MULTIPLIER 

SPEC STEP MULTIPLIER 

OBJECTIVE ^ 

MAGNIFlCATinN 

If SOLUTION 

Q 

Q 9_ 

2 2 

__1 0 

0 

2 

J Q Q_ 


cnne 


FILM NUMBER 


. ... MONTH 

DAY 

V r AP 

n 

1 

1 

1 

1 

1 

0 

t 

1 * 

7 A 


CODE 


CAL 

STRIP NUMBER 

■ 

STEP tow PENS ITT 

STEP NIGH DENSITY 


X 

__JL_ 

5 

3_ 

Q 

5 

0 

l 

i i 

<l, a. 




TABLE II.- STEP WEDGE DENSITY VERSUS "BCD" NUMBER TABLE 


Density 

BCD No. 

0.05 

71 

.20 

85 

.36 

97 

• 50 

111 

.6k 

12ir 

.79 

138 

.9^ 

151 

1.09 

163 

1.23 

1 

1.38 

lh 

1-53 
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^CD numbers go up to 175 and then start over again at 1, 


ro 



















Center trace 


Starting point 
Border trace 




Fig. 1. - Sample photograph and Den plots 
of line 51 and bottom border 
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Fig. 2. - Pictorial representation of camera system angles 



(a) The densities from each 
trace are plotted from right 
to left. Trace one is at the 
bottom. The computer picture j 
is thus an upside down version 1 
of the photograph shown in ' 






